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We present extensive results and analysis of energy and angular distributions of diffuse UHE u e , 
Vp, and v T fluxes propagated through earth, with and without augmentation of the standard model 
interactions by low scale gravity. With propagated fluxes in hand we estimate event rates in a lkm 3 
detector in ice with characteristics of ICECUBE. We determine that, at 0.5PeV energy threshold, 
there is a significant difference in the ratios of down shower events to upward muon events between 
the standard model and the low scale gravity cases with ITeV and 2TeV mass scales. The same 
is true for energy threshold at 5PeV. Though the difference is large in all flux models, statistical 
significance of this difference depends on the flux models, especially at 5PeV and above. Both 
flavor assumptions, v B , f M , v T :: f, 2, and v e , z/ M , v T :: 1, 1, 1, and all flux models show large 
differences. Though rates of tagged events are low, we find that v T regeneration by r decay may 
play an important role in disclosing deviations from standard model predictions at energies in the 
neighborhood of f PeV for ITeV-scale gravity, for example. We emphasize those analyses whose 
sensitivity to new physics is independent of the flux model assumed. 

PACS numbers: 96.40.Tv, 04.50.th, 13.15.+g, 14.60.Pq 



I. INTRODUCTION 

The pursuit of high and ultra-high energy neutrinos 
has greatly intensified over the past decade as more 
and more neutrino telescopes have entered the search. 
Though the observation of MeV-neutrinos emitted from 
SN 1987a is over 15 years old 0,0, there is still no firm 
candidate for TeV, PeV or EeV neutrinos of galactic or 
extra-galactic origin. Yet there is good reason to expect 
a neutrino flux exists in this energy regime because of 
the great success of air-shower detectors in building a 
detailed record of cosmic rays with these very-high to 
ultra-high energies 0. The photons or nuclear parti- 
cles that are generally believed to initiate the observed 
shower are accompanied by neutrinos with similar ener- 
gies, in most models of the high-energy particle emission 
by the sources. In any case, the neutrinos emitted by 
production and decay of pions by the highest energy pri- 
mary cosmic rays as they interact with the cosmic mi- 
crowave background, the so-called GZK |j] neutrinos, 
should be present at some level at ultra-high energies 
|{|, regardless of the mechanism responsible for produc- 
ing the observed cosmic rays H, 0, Hi ■ Even if there are 
no super- GZK neutrinos, there are a number of models 
HEmEHQEHEi, which Predict the exis- 
tence of neutrinos in the PeV-EeV range. By choosing 
several contrasting flux models and usi ng e nhanced cross 
sections from low scale gravity 0, 0, 0| , we look for 
new physics effects that are relatively independent of flux 
models. 

The expanding experimental capabilities and the 
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strong theoretical interest in understanding the physics 
of astrophysical sources and particle interactions of the 
highest energy cosmic rays makes it imperative to study 
all aspects of the neutrino observation process. The num- 
ber of groups reporting limits on fluxes and projecting 
improved limits with expanded data sets or with new fa- 
cilities is impressive. In the range 1 TeV to 1 PeV, the 
AMANDA pjjj, Frejus Hj, MACRO and Baikal H| 
experiments have reported limits on neutrinos from as- 
trophysical (non-atmospheric) sources. In the range 1 
PeV to 1 EeV, AGASA [H, AMANDA [H, Fly's Eye 
[H, and RICE have all reported limits. Above 1 
EeV, AGASA, Fly's Eye, GLUE H| and RICE all put 
limits on the flux that extend up into the GZK range. 
The upper limits are getting interestingly close to the 
predictions of several models and actually below the pre- 
dictions in several cases. The situation is heating up and 
will get hotter as the experiments like AUGER, which is 
already reporting preliminary results on air showers p9| 
and ICECUBE |3(j] are fully operational. Meanwhile, ex- 
panded data sets and improvements in sensitivity in ex- 
periments like RICE will continue to search and to push 
down on limits until the first UHE neutrinos are observed 

m 

These detection capabilities that have been achieved 
and will be improved and expanded in the next few years 
have direct impact on particle physics. The detection 
estimates, upon which limits are based, all rely on the 
extrapolation of neutrino cross sections well beyond the 
currently measured energy r ange . Is QCD correctly pre- 
dicting these cross sections 32]? Is there new physics 
that enhances neutrino cross sections at high energies 
[33? What is the effect of new neutrino interactions 

HSIHIillillHli^or neutrino mixing [UIHIII 
on the expected rates of detection in various telescopes? 
Clearly there is ample motivation for examining the con- 
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sequences of various combinations of assumptions about 
the physics governing the cross sections and the assump- 
tions about the flavor composition of the astrophysical 
flux of neutrinos. What, if any, are the observable dis- 
tinctions among the various possibilities of flux and in- 
teraction characteristics? These questions and the exper- 
imental prospects for answers motivate this work. 

There is considerable published work on r-neutrino 
propagation through earth in standard model fSM) us- 
ing analytic and computational tools [43l lil . lia . Ea . ItH 
lid I4H [5fl l5lj| in the scenarios v e , v^, v T :: 1, 2, and 
j/ e , v^, v T :: 1, 1, 1, and some analytical and computa- 
tional work on neutrino propagation in low scale gravity 
(LSG) models has also been done in the v ei , v T :: 1, 
2, scenario |34L l52| . A detailed study has not been 
done in the v e , v^, v r :: 1, 1, 1 scenario in LSG mod- 
els. In this paper we solve, using Runge-Kutta method 
[53|. the coupled differential equations for the four lep- 
tons v e , v^, v T , and r in both of the above scenarios, in 
SM and LSG models. For cross section calculations, we 
use Gaussian and monte carlo integration methods |53l ] 
with CTEQ6-DIS parton distributions jH^I- Our results 
confirm significant regeneration effect due to taus in the 
SM as already shown by several authors 0, 01 • How- 
ever, as we will see the regeneration due to taus is not as 
significant in LSG models. Also, by comparing results of 
[4a. It^ | , one finds that electromagnetic (EM) losses of r 
are not making a significant difference in the SM fluxes 
of v r around IPeV, hence, we do not include EM losses 
in our work here. As we will see in the next section, EM 
losses are not important at all in LSG models. 

In Section 2 we talk about cross sections and interac- 
tion lengths in SM and LSG ; Section 3 gives the for- 
malism for neutrino propagation through the earth; in 
Section 4 we show our results for different neutrino mod- 
els and discuss them; in Section 5 we develop formalism 
for event rates calculation and in Section 6 we show and 
discuss our results for event rates; Section 7 gives the 
summary of our results and the conclusion. 



we will use NC and CC cross sections for all of them. 
The differences between r and v weak cross sections are 
not important because r decay is the only dominant pro- 
cess for t's up to energies 10 8 GeV, and at higher ener- 
gies the differences are not significant anyway (55|. A 
detailed discussion of SM z/ M — nucleon cross sections is 
given in Ref.[56||. We have used CTEQ6-DIS parton dis- 
tributions. Details of calculation for LSG model cross 
sections are given in Ref.. IMOS, where CTEQ4-DIS 
parton distributions [5^ were used. The difference be- 
tween CTEQ4-DIS and CTEQ6-DIS cross section calcu- 
lations is not significant. We outline the LSG calculation 
here so that the presentation is reasonably self contained. 

The classical gravity Schwarzschild radius rs(*/s) is 
the dominant physical scale when the collision energy is 
large compared to the Plank mass. At impact parameters 
smaller than rs , we use the parton-level geometrical cross 
section 591 
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For values of the classical impact parameter, b, larger 
than r s we use the contributions to the amplitude in the 
eikonal approximation. In Eq. ^ rs is the Schwarzchild 
radius of a 4 + n dimensional black hole of mass Mbh = 

vs m, 
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where vs is the neutrino- parton CM. energy, and M is 
the 4+n-dimensional scale of quantum gravity. Multiply- 
ing by the parton distribution functions, fi(x,q), choos- 
ing q at a value characteristic of black hole production 
and integrating over momentum fraction x, gives the es- 
timate 



dxa 



■BH{xs)j % {x,q). (3) 



II. CROSS SECTIONS AND INTERACTION 
LENGTHS 

Before evaluating the equations of propagation for the 
four leptons v e , v^, v Tl and r through the earth in 
next section, we need to calculate their cross sections on 
isoscalar nucleons (N = ^rp), where n stands for neu- 
tron and p for proton. We need to calculate their neutral 
current (NC) and charged current (CC) weak interaction 
cross sections in SM and eikonal (EK) and black hole 
(BH) cross sections in LSG models. The LSG models do 
not discriminate among different particles; they are the 
same for all the four leptons v ei v^, v T , and r. The SM 
total cross sections for different flavors are also the same 
within a few percent at the ultra high energies we are 
interested in here. We will assume the SM cross sections 
are the same for all the four leptons v e , v^, v T , and r, and 



We take x m i n = M 2 /s and q = ^Js. The dependence of 
o~vN^bh(s) on the choice of Xm in and the treatment of 
q is discussed in Ref. [6l|, [63, 163, [64J ■ 

For the input amplitude to the eikonal approximation, 
referred to as the Born amplitude, we choose 

iMsom = 2^ "772" "IT 2 ' 4 

j 3 

where c is the gravitational coupling strength, c = 
(M/Mp) 2 and M P = 2.4 x 10 18 GeV is the reduced, four 
dimensional Planck mass. Here q = t is the usual 
lepton momentum transfer. The index j must include 
the mass degeneracy for the jth K-K mode mass value. 
The sum, which can be well approximated by an integral, 
must be cut off at some scale, generally taken to be of 
the order of M. The transverse Fourier transform of the 



Born amplitude produces the eikonal phase as a function 
of impact parameter b, 

i f d 2 q 

% ( s ' & ) = 2^ / 4^2 ex P(*q- h)iM B om- (5) 

Evaluating the integral over q and representing the sum 
in the Born term by an integral, one finds the ultraviolet- 
finite result 
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The eikonal amplitude is then given in terms of the 
eikonal phase by 

M = —2is J d 2 6exp(iq • b) [exp(ix) — 1] 

= -Utts J dbbJo(qb) [exp(i X ) - 1] . (8) 

The eikonal amplitude can be obtained analytically |65l 
0,113 m the strong coupling qb c >> 1 and weak coupling 
qb c << 1 limits. 

For strong coupling, the stationary phase approxima- 
tion is valid, yielding 
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and b s — b c (qb c /n)~ 1 '( n+1 \ In the cross section calcula- 
tion, we set the amplitude equal to its value at q = l/b c 
for values of q that are less than l/b c , since the small q 
region makes negligible contribution to the cross section. 

We assume that the black hole cross section is the dom- 
inant one for q > 1/rg. The eikonal cross section is cut 
off at this value of q, since it is not expected to be reliable 
for values of q larger than 1/rs in any case. 

In Fig.l. we plot the SM and LSG (models with mass 
scale ITeV and 2TeV; n = 6) neutrino-isoscalar nucleon 
cross sections. For SM we plot neutral current (NC), 
charged current (CC), and total ( NC + CC) cross sec- 
tions. For LSG we plot eikonal (EK), black hole (BH), 
and total (EK+BH) cross sections. We see obh is larger 




FIG. 1: ^-isoscalar nucleon cross sections (cm 2 ) vs energy 
(GeV): for low scale gravity (LSG) models, with number 
of extra dimensions n=6, we plot eikonal (EK, dashed line), 
black hole (BH, dotted line), and total (EK+BH, solid line) 
cross sections; for standard model (SM) we plot neutral cur- 
rent (NC, dashed line), charged current (CC, dotted line), 
and total (NC+CC, solid line) cross sections. 




FIG. 2: Interaction lengths (km) vs energy (GeV): interac- 
tion length dnt = a (E)N AP ' w here a(E) are plotted in Fig. 1, 
Na is Avogadro's number, and p is the material density; we 
choose p = 8g.cm~ s to make some comparisons with Ref. |4fi| | 
. For low scale gravity (LSG) models, with number of extra 
dimensions n=6, we plot eikonal (EK, dashed line), black hole 
(BH, dotted line), and total (EK+BH, solid line) interaction 
lengths. For standard model (SM) we plot neutral current 
(NC, dashed line), charged current (CC, dotted line), and 
total (NC+CC, solid line) interaction lengths. We also plot 
r-decay length (thick solid line). 



then oek for our case (n — 6). However, in n = 3 case 
not shown here, the reverse is true. Our results turned 
out to be only marginally sensitive to the number of di- 
mensions, so we choose to work with n — 6, for which 
the bound on the scale M is the weakest. 

The interaction length in a material with density p is 
defined here as Ci n t — ( N Apa )i where Na is Avogadro's 
number and a is the cross section for the interaction. 
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Fig. 2. gives the interaction lengths in SM and LSG as 
well the r-decay length. We set p — 8.0gcm~ 3 to make 
some comparisons between the interaction lengths and 
the electro mag netic (EM) ranges of taus and muons given 
in the Ref . |45| . using this value of p. 

Which particles do we need to include in our propaga- 
tion of neutrinos through the earth? We have six can- 
didates which might be coupled with each other: v ei v tll 
v T , and their three leptonic partners. One can exclude 
electrons from this list, both in SM and LSG, because 
they shower electromagnetically before they produce any 
v e through the CC interaction at energies of our interest 
here. The case for p's and r's needs some attention be- 
cause it is different in SM and LSG. In SM we can ignore 
p's but not r's because of the comparatively smaller de- 
cay length and larger EM ranges of taus. For LSG, muons 
play some role in the propagation through the earth at 
high enough energy as explained below. The effect is 
quite small, however. 

If we look at the interaction lengths and r-decay length 
in Fig.2, we find that: (i) the LSG1 model (LSG model 
with mass scale ITeV) interaction lengths become smaller 
than the r-decay length for E > 10 8 GeV. This implies 
the regeneration effect due to taus will be suppressed in 
LSG models (ii) if we look at the EM ranges of taus 0] , 
we find that the LSG1 model interaction lengths become 
smaller than the tau EM range for E > 10 s GeV. This 
gives a reason for not including, in LSG, the EM energy 
losses of taus in our propagation of neutrinos, coupled 



with taus, through the earth, (iii) If we look at the EM 
ranges of muons |4^, we find that even the muon EM 
energy losses are not important for E > 10 s GeV. (iv) 
Given the above reasons, interestingly, taus and muons 
become almost identical in LSG for E > 10 8 GeV. This 
means one may have to treat muons and taus on equal 
footing in the propagation of neutrinos through the earth 
for E > 10 8 GeV in LSG1 model and for E > 10 9 GeV in 
LSG2 model (LSG model with mass scale 2TeV). How- 
ever, in the present work, we do not include muons in the 
propagation equations because we are looking at neutri- 
nos around IPeV here, and we expect the coupling of 
muons with the via CC interaction for E > 10 8 Gey 
will not affect the neutrino flux much around IPeV by 
feed down; the muon decay length, being so large in con- 
trast to that of r, will play no role in regeneration of 

mi- 



ni. NEUTRINO PROPAGATION THROUGH 
EARTH 

Here we discuss the coupled propagation of v e , v^, v Tl 
and r through earth. We do not include the EM energy 
losses of taus for the reasons discussed in t he p revious sec- 
tion. Suppose we have a differential flux 83] F l (E,x,9) 
of lepton of species i at the surface of earth, then the 
transport equation for each of the four leptons is 



dF l (E,x, 
dx 
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where the first two terms give the loss and the last two 
terms give gain of the flux per unit length in the same 
energy bin E. Avogadro's number Na times the density 
p(x, 8), gives the number of target nucleons per unit vol- 
ume at the nadir angle 9 and distance x in the earth (se e 
Fig. 3.). We use the earth density model from Ref. |68| . 
a\ (E) is the total cross section for a lepton of flavor i to 
interact with a nucleon and be expelled from the energy 
bin E: 



4(E) = ahdE) + a* NC (E) + a* BH (E) + a* EK (E), 

where we use the same (j\(E) for all the four leptons v e , 
v^, v Tl and r (LSG cross sections are the same for them; 
see Section 2 for discussion on SM cross sections.). The 
second term in Eq^| gives the loss due to decays. It is 



zero for the neutrinos, and for taus 

Q ec (E) = 7 cT, (13) 

where 7 = — is the Lorentz factor, T is the mean life 
time of taus and c is the speed of light in vacuum. The 
third term in Eo 1 121 gives the gain in the flux of species i 
in the bin E, resulting from interaction of the species j 
at E > E, and 

^ da^jE'.E) = da^(E',E) da 3 ^^ (E , E) 

^ dE dE dE 

j 

dE 

In the above equation, there is no need for the second 
term on the right hand side for v e and propagation 




FIG. 3: Drawing of the earth: 9 is the nadir angle and Re is the radius of the earth. Arrows outside the earth represent the 
down flux. Upflux is the flux coming through the earth. 



equations as we are not keeping track of electrons and 
muons for the reasons given in the last section. How- 
ever, we keep this term for the v T and r equations as 
we are propagating taus along with the neutrinos. For 
the reasons given earlier, we use the same NC, CC, EK 
differential cross sections for all the four leptons in the 
above equation. The fourth term in Eq^|is the gain in 
flux of species i in the bin E due to tau decays. This term 
is zero in the tau flux equation. For v e , v^, and v T we 
consider the corresponding decay channels of taus as for- 
mulated in [46l 2000] :(i) t — > v r \w^ (ii) t — > u T eu e , (iii) 
t — * u t tt, (iv) t — > v T p, (v) t — > v r al, and (vi) r — > v T X . 
These decays have branching ratios of 0.18, 0.18, 0.12, 
0.26, 0.13, 0.13 [6|j, respectively. All the decays give a 
v T . The first two decays couple the tau propagation with 
v e and propagation. The last decay includes the rest 
of the hadronic decays not specified in (iii) through (vi). 
For a review, see |70tl7l| . 

IV. RESULTS AND DISCUSSION FOR FLUXES 

Below we show results for total fluxes, including all 
the neutrino species and their respective antineutrinos, 
unless defined otherwise. We plot upward fluxes instead 
of up-to-down flux ratios to find the region in the (E, 8)- 
space to compare the SM and LSG models in terms of ab- 
solute flux differences. Larger fluxes mean more events, 
and larger difference in the number of neutrinos in SM 
and LSG means better chances to differentiate between 
the models. The ratio plots are not always helpful for 
that purpose because they do not show us the actual 
number of neutrinos and the flux difference of the SM 
and LSG. We present plots for AFi(E, 9), AF 2 (E : 9) in 
the (E, 6>)-space, where, 

A-Fi (E, 9) = F SM (E,9) - F LSG1 {E,B), (14) 

and, 

AF 2 {E,6) = F SM (E,e) - F LSG2 {E,6). (15) 

Here F L sgi(E,9), F LS g2{E,9) are the total upward 
fluxes in the low scale gravity models with mass scale 
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FIG. 4: Upward neutrino flux (fem _2 .j/r _1 .sr _1 )integrated 
over energy E > 0.5PeV vs nadir angle 6 (degrees) for R. J. 
Protheroe input flux model: SM upward neutrino flux Fsm{9) 
, and flux differences AFi(8), and AF2(6), as defined in Eqs. 
lf4l and lT51 are plotted for two scenarios: v & : Vy, : v T :: 1 : 1 : 1 
(solid lines), and v e : : v T :: 1 : 2 : (dotted lines). 



ITeV and 2TeV respectively, and number of extra dimen- 
sions n = 6. Fsm{E,9) is the upward flux in standard 
model. AF\{E), AF 2 {E), AFx(9), and AF 2 {8) are de- 
fined in the same manner as AFi(E, 8) and AF 2 {E,9) 
are defined above. 

Below we show our analysis for the neutrino flux mod- 
els due to R. J. Protheroe [nl, K. Mannheim (B) [l^ . 
Waxman Bahcall E3, SDSS 10, and 1/E generic model. 
Though we did not show it here, we also looked at at- 
mospheric and galactic neutrinos [T^ around 0.5PeV . 
However, the up fluxes in this case are ignorable as com- 
pared to the above extragalactic flux models: the galactic 
up flux is more than an order of magnitude smaller than 
the up flux from any of the flux models considered here; 
the atmospheric up flux is more than two orders of mag- 
nitude lower than even the galactic up flux. The reason is 
quite simple: one should expect larger feed down in the 
upward model fluxes considered here; these extragalac- 
tic models have much larger fluxes above IPeV than the 
atmospheric or galactic fluxes. As a result, the model 
fluxes at higher energies give larger feed down at 0.5PeV 
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FIG. 5: Upward v T flux (fcm _2 .j/r~ 1 .sr~ 1 ) integrated over 
energy E > O.bPeV vs nadir angle 8 (degrees) for R. J. 
Protheroe input flux model: SM upward v T flux Fg T M (6) (solid 
line), and v T flux differences AF" T (6), and AF^lff) (dot- 
ted lines), as defined in Eqs. 1141 and Il5l for the scenarios: 
v e : v,j, : v T :: 1 : 1 : 1. 
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FIG. 6: Same as Fig. 5, but with K. Mannheim (B) flux 
model. 



while the atmospheric and galactic up fluxes are so small 
there that they give essentially no feed down. This is es- 
pecially true in SM because of the strong feed down effect 
per incident v T due to tau decays, which may make the 
extragalactic model up fluxes, in SM, large enough to be 
detectable even at 0° nadir angles around 0.5PeV . 

We will first show some plots (Figs. 4-14) to explore 
(E, 0)-space of the neutrino fluxes, and finally we will 
give two tables of numbers for different neutrino flavor 
flux ratios. Three types of plots are shown below: 

(i) Fluxes integrated over energy vs nadir an gle: 
Figs. 4 and 5 refer to our example of R. J. Protheroe |ll| 
model. Fig. 4 shows the angular distribution of the to- 




FIG. 7: Same as Fig. 5, but with Waxman Bahcall (WB) flux 
model. 
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FIG. 8: Same as Fig. 5, but with SDSS flux model. 



tal SM flux Fsm(0) integrated over energy E > 0.5PeV, 
along with the flux differences AFi(E) and Ai^i?) as 
defined in the Eqs^J an( l 1^1 The qualitative features 
of the Figs. 4 and 5 are the same in every flux model, so 
we do not show figures for the other flux models corre- 
sponding to Fig. 4. However, we show the corresponding 
figures for Fig. 5 for the four models (R. J. Protheroe 
Ml, K. Mannheim (B) [lp, Waxman Bahcall P3, SDSS 
|9|) because our emphasis here is on r-effects, and we 
wish to show their qualitative features are independent 
of the flux model. As mentioned earlier, the atmospheric 
and galactic ff^ up fluxes are not significant as compared 
to the four extragalactic source models at these energies, 
so we do not consider them here. 

We plot the upward fluxes in two scenarios for the ini- 
tial flux: v e , v^, v T :: 1, 1, 1 (solid lines) and v e , v^, 
v T :: 1, 2, (dotted lines) . The two scenarios cor- 
respond to Va — > v T oscillations and no oscillations in 
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space, respectively. The following observations are com- 
mon to all models and are relevant to Fig. 4: (i) the 
difference is largest in number around 75°-85°; however, 
as we will see in 3D plots in (E, f?)-space, the maximum 
of the difference shifts to lower angles at lower energies, 
(ii) around 30° nadir angle, one can clearly see the effect 
of the core — stronger suppression, (iii) the difference in 
flux between SM and LSG models is larger at any angle 
in the v e , v^, v T :: 1, 1, 1 scenario than in the no tau 
scenario. This is expected because of the stronger tau 
regeneration effect in the SM as compared to LSG; in 
LSG, the black hole cross section, being the largest of all 
as shown in Fig. 1, suppresses regeneration due to any 
process. This observation leads us to concentrate on tau 
fluxes only, as discussed below. 

In Figs. 5-8 we show the same plots for v T as was shown 
for the total fluxes. We plot these figures for all the four 
models: R. J. Protheroe [lj, K. Mannheim (B) pf . 
Waxman Bahcall E3, and SDSS @. In these fi gures, 
the solid line is the v r flux in SM and the two dotted 
lines are the difference fluxes (see Eqs^J ^3 ) for v T . 
We clearly see that in LSG1 model, regardless of the 
flux model, difference in v T fluxes is more than 50 % 
of the total difference due to all neutrino species. For 
example, the v T difference AFi [E) (Fig. 5) for Protheroe 
model at 0° nadir angle is around 3000, while the total 
difference is around 4000 (Fig. 4), km -2 yr -1 sr -1 . This 
behavior is independent of the flux models — about 3/4 
of the total difference is due to v T only. This may be 
useful as ICECUBE is expected to differentiate between 
neutrino flavors around 500 TeV |7^| . 

(ii) Flux integrated over nadir angle vs energy: 

The fluxes integrated over nadir angle do not show as 
much detail as the ones integrated over energy. This is 
because the integrated flux gets dominant contribution 
from nadir angles around 90° where the chord length of 
the earth is not long enough to make the difference be- 
tween SM and LSG prominent. For this reason, it suf- 
fices to show plots for only one flux model , chosen to 
be the R. J. Protheroe model [llj; also the qualitative 
features of these plots, like the earlier plots, are model 
independent so it is not important to show the figures 
for all the models. We see in Fig. 9: (a) Even at en- 
ergies E — 0.5PeV, the LSG1 and SM flux models are 
distinguishable; (b) in the v e , v^, v T :: 1, 1, 1 scenario, 
AF\{E) and AF2(E) are larger at lower energies, how- 
ever, in contrast in the v e , v^, v r :: 1, 2, scenario, 
they decrease with decreasing energy. This contrasting 
behavior of the two scenarios begins around IPeV for 
LSG1 and around 5PeV for LSG2. This is expected be- 
cause the stronger feed down effect due to taus causes 
Fsm(E) in the z/ e , v T :: 1, 1, 1 scenario to increase 
faster with decreasing energy than in the other scenario; 
however, Flsgi(E) and Flsg2(E) are not as different in 
the two scenarios because they are not as sensitive to tau 
regeneration. Hence, the over all result is decreasing flux 
difference for the v e , v^, v T :: 1, 2, scenario and increas- 
ing flux difference for the other scenario. This happens at 
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BTeV) 

FIG. 9: Upward neutrino flux (km~ 2 .yr~ x .TeV -1 ) integrated 
over nadir angle vs energy E (PeV) for R. J. Protheroe in- 
put flux model: SM upward neutrino flux Fsm{E) , and flux 
differences AFi(E), and AF2(E), as defined in Eqs. 1141 and 
1151 are plotted for two scenarios: v e : Vp : v T :: 1 : 1 : 1 (solid 
lines), and v K : : v T :: 1 : 2 : (dotted lines). Also shown 
is the downward R. J. Protheroe model flux, integrated over 
nadir angle (dashed line). 

a lower energy in the LSG1 model due to the lower energy 
scale for LSGI (Figs.l, 2); (c) another important obser- 
vation is that AFi(E) is almost equal to Fsm(E) after 
3PeV, which means if we want to differentiate the two 
models on the basis of the nadir angle integrated event 
rates, the best region in energy may be around 3PeV, 
if the detector has large enough efficiency to detect this 
flux. However, obtaining larger fluxes for better statis- 
tics requires looking at lower energies. Though the model 
fluxes around energies as low as 100 TeV are larger, the 
percent difference between SM and LSG fluxes becomes 
smaller and smaller at energies below 0.5PeV where SM 
cross sections are dominant; the total upward flux rises 
much faster than the flux difference, making it hard to 
differentiate between the two models. The atmospheric 
background is also larger at these energies, hence we did 
not find it interesting to show the fluxes below 0.5PeV. 

In Fig. 10 we plot the nadir angle integrated flux of v T 
only. If we compare Fig. 9 and 10, we come up with the 
similar answer as we did for the flux integrated over en- 
ergy (AF" T (E)): the upward v T flux difference, AF" T (6), 
around 0.5PeV is almost 3/4 of the total upward flux dif- 
ference AFi (6) . This again gives one hope that the sig- 
nals of low scale gravity may appear even around 0.5PeV. 

(iii) Plots of flux as a function of both energy 
and nadir angle: Figs. 11-14 give the complete detail of 
the fluxes in the (E, #)-space for our flux example of R. J. 
Protheroe [llj. In Figs. 11 and 12 we plot the total up- 
ward flux in SM , Fsm{E,9), and the total upward flux 
difference AFi(E,6), respectively; Figs. 13 and 14 have 
similar plots for the v T . We can see in these plots: (a) The 
difference is the largest around 80°, however, it is still in- 
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Integrated over nadir angle- 



AF ',(8) 



t — 1 — i — 1 — i — 1 — r 



PR (V : V :V : : 1 : 1 : 1 scenario) 




FIG. 10: Upward v T flux (fcm _2 .j/r~ 1 .Tel /_1 ) integrated over 
nadir angle vs energy E (PeV) for R. J. Protheroe input flux 
model: SM upward v r flux F$ T M (E) (solid line), and v T flux 
differences AF" T (E), and AF^(E) (dotted lines), as defined 
in Eqs. 1141 and [T51 for the scenarios: v e : : v T :: 1 : 1 : 1. 
Also shown is the downward v T R. J. Protheroe model flux, 
integrated over nadir angle (dashed line). 



creasing even at 0.5PeV. Again, this may be surprising 
at first glance because below IPeV there is no significant 
contribution to the cross sections from LSG. However, 
the reason is simply that the cross section at a given 
energy will affect the neutrino flux at equal and lower 
energies due to feed down. Keeping this in mind, we 
can argue that the flux around IPeV or below gets more 
feed down from higher energies in SM because LSG black 
hole cross section, being the largest of all the cross sec- 
tions, suppresses the feed down effect due to any process, 
(b) At higher energies the flux difference peaks at higher 
angles e.g. around 80°; the peak shifts to lower angles 
at lower energies. One may argue that the peak should 
always occur at the lowest nadir angle because the neu- 
trinos will have more interactions as they pass through 
earth with larger chord lengths, and hence the SM and 
LSG models' interaction will cause the flux differences 
to become larger and larger at lower angles and higher 
energies. However, this does not happen because the in- 
put fluxes at higher energies are so small that at lower 
nadir angles all the flux is either absorbed or fed down 
to lower energies; that is why we see the flux difference 
peak shifts towards 0° nadir angle at lower energies: feed 
down effect makes the difference, between SM and LSG, 
at higher energies appear at lower energies, (c) The v r 
plots in Figs. 13, 14 show us that the major contributor 
of the difference between SM and LSG is the v T at lower 
energies. It contributes almost 3/4 of the total difference 
around 0.5PeV and around 1/3 at lOPeV. This feature is 
best seen in these full (E, #)-space plots. Around energies 
0.5PeV, v T plays important role in probing new physics. 
At energies around lOPeV and higher, v T behaves more 
like v & and v^; this is because around these energies the 
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SM upward neutrino flux Fsm(E,9) 
(km~ 2 .yr^ 1 .TeV~ x .sr _1 ) vs energy E (PeV) and 
nadir angle 9(deg), for R. J. Protheroe input flux model in 
the scenario: v e : u u '■ v T :: 1 : 1 : 1. 



PR (v :v : v : : 1 : 1 : 1 scenario) 
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Upward neutrino flux difference AFi(E,9) 
(km~~ 2 .yr~ Y .TeV~ x .sr" 1 ) vs energy E (PeV) and nadir 
angle 8(deg), for R. J. Protheroe input flux model in the sce- 
nario: i>» : v a : v T :: 1 : 1 : 1. 



feed down due to taus, from even higher energies, is not a 
big effect both in SM and LSG. This is a result of larger 
tau decay lengths, smaller interaction lengths (see Figs. 
1, 2), and smaller fluxes at higher energies. 

Summarizing, we see that the above analysis discloses 
the flux structure in(E, #)-space: (i) Around 0.5PeV, the 
flux difference peaks in the 40 — 60° region. The larger 
angles tend to wash out the difference between LSG and 
SM. (ii) Around 5PeV, the difference peaks in the 75-80° 
region, (iii) At higher energies, one will have to look at 
even larger nadir angles to get any detectable up flux. 
In Tables HI and ITT1 we give different flux ratios Rl, R2, 
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PR (v -v :v : : 1 : 1 : 1 scenario) 




FIG. 



13: SM upward v T flux F S m(E,6) 

(km~~ 2 .yr^ 1 .TeV~ x .sr" 1 ) vs energy E (PeV) and 
nadir angle #(deg), for R. J. Protheroe input flux model in 
the scenario: i/ e : : v T :: 1 : 1 : 1. 

PR (V :v :v : : 1 : 1 : 1 scenario) 




FIG. 14: Upward v T flux difference AF^{E,8) 
{kmT 2 .yr^ 1 .TeV~ x .sr _1 ) vs energy E (PeV) and 
nadir angle #(deg), for R. J. Protheroe input flux model in 
the scenario: v e : v a : v T :: 1 : 1 : 1. 



and i?3 defines as, 

R1 _ F total M R2= Fu T (np) 



R3 



in the v e , v, 



F total (down) 

*V>P) 

F, e (up)+F„ (up)' 



F„ T (down) ' 



(16) 



: 1, 1, 1 scenario, for SM, LSG with 
M=2TeV (G2), and LSG with M=lTeV (Gl). These 
ratios reveal some interesting and useful features: (i) If 
we compare R2 and i?3 with Rl in the Tables d and HH 
we find that v T flux indeed behaves differently in SM 
and LSG. (ii) In LSG, at higher energies, v e , v^, and v T 
tend to become identical as expected from larger decay 



TABLE I: Flux ratios of the fluxes integrated over energy 
E > 0.5PeV at a fixed nadir angle 9 = 45°, for SM, LSG 
with M = 2TeV (G2), and LSG with M=lTeV (Gl); Rl, R2, 
and R3 are defined in Eqs. 1161 



MODEL 



WB 
M(B) 
PR 
SDSS 



Rl 



SM 
0.18 
0.15 
0.17 
0.18 
0.10 



G2 

0.17 

0.12 

0.16 

0.18 

0.05 



Gl 

0.12 

0.06 

0.09 

0.12 

0.02 



R2 



SM 
0.25 
0.29 
0.27 
0.24 
0.22 



G2 

0.23 

0.20 

0.24 

0.24 

0.09 



Gl 

0.14 

0.07 

0.12 

0.15 

0.02 



RS 



SM 
0.85 
1.61 
1.13 
0.83 
2.84 



G2 

0.82 

1.26 

1.03 

0.82 

1.60 



Gl 

0.66 

0.73 

0.73 

0.66 

0.80 



TABLE II: Flux ratios of the fluxes integrated over nadir angle 
6 at a fixed energy E = 5PeV , for SM, LSG with M = 2TeV 
(columns G2), and LSG with M = ITeV (columns Gl); Rl, 
R2, and f?3 are defined in Eqs. 1161 ) 



MODEL 


Rl 


R2 


R3 




SM 


G2 


Gl 


SM 


G2 


Gl 


SM 


G2 


Gl 


WB 


0.20 


0.18 


0.02 


0.25 


0.22 


0.02 


0.71 


0.65 


0.53 


M(B) 


0.28 


0.22 


0.02 


0.45 


0.30 


0.02 


1.12 


0.84 


0.56 


PR 


0.21 


0.19 


0.02 


0.29 


0.23 


0.02 


0.80 


0.69 


0.54 


SDSS 


0.19 


0.18 


0.02 


0.23 


0.20 


0.02 


0.66 


0.62 


0.53 


1/E 


0.43 


0.24 


0.02 


0.82 


0.35 


0.02 


1.72 


0.93 


0.57 



length of taus at higher energies and larger black hole 
cross sections hence smaller feed down of v T from tau de- 
cays; this effect is indeed seen for LSG (ITeV) (see Gl 
in Table [nj; Rl and R2 tend to become equal and R3 
tends to become 0.5 as expected (see Gl in Table lTT|l . The 
same will be true for LSG (2TeV) at even higher ener- 
gies, (hi) No matter which neutrino flux model is correct, 
even at energies as low as 0.5PeV, we see a clear differ- 
ence between SM and LSG (ITeV) model (Table HJ; at 
higher energies, E > lOPeV, there should be a difference 
between LSG(2TeV) and SM based on their their com- 
parison at 5PeV (Table nj). Our data at lOPeV which 
is not shown here, implies that, based on the flux ra- 
tios, the % difference between LSG(2TeV) and SM at 
lOPeV is bigger than the one between LSG (ITeV) and 
SM at 5PeV.(iv) Isolating v T should help differentiate 
between SM and LSG, and may help to differentiate be- 
tween different neutrino flux models given SM dynamics. 
For example, in Tables [I] and [HI though Rl is similar 
for different flux models in the SM case, i?3 shows some 
significant variation. 

However, results of propagated WB and SDSS models 
are close to each other; they are similar because, though 
WB down flux at lower energies is much smaller than 
SDSS, WB flux falls like E~ 2 while SDSS goes like E~ 3 
at higher energies giving a stronger feed down effect for 
the former, hence, their up-down flux ratios tend to be 
the same; the big difference between these models is that 
WB has much weaker flux at energies 0.1-10PeV (Figs. 6 
and 7 ) . (v) The difference among flux models are largely 
washed out by the LSG dynamics at higher energies. For 
example, in Table [H] Rl and R2 are the same for all the 
flux models given here when the LSG scale is ITeV. As 
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mentioned above, to probe higher LSG scales, one must 
go to higher energy data which shows some sensitivity 
upto M = 5TeV [a|. 



V. EVENT RATES 

Next we outline the formalism to calculate event rates 
for showers, muons, and taus. Our formalism adds some 
refinement to the event rate estimates |42l Fij . In pre- 
senting event rates, we take our theoretical "ICECUBE- 
like" detector to be 1 km 3 of strings of optical modules 
deployed with 125m horizontal spacing. Because rates 
of events depend on flux and cross section, the extra de- 
pletion of upward flux in the LSG models is somewhat 
compensated by the increased interaction probability of 
each neutrino that penetrates the detector effective vol- 
ume. 

(i) Shower rates: Neutrino-nucleon interactions at 
PeV energies and above initiate electromagnetic and 
hadronic showers which may produce a detectable ra- 
dio or optical Cherenkov signal in detectors like RICE, 
AMANDA, ICECUBE. For shower rates in SM, we in- 
clude both CC and NC interactions of v e and v T but only 
NC interactions for (CC interaction in this case gives 
muons which can be detected directly in detectors like 
AMANDA and ICECUBE). For LSG, we include both 
eikonal and black hole cross sections. For shower rates 

Rshower h*l LSG, 

oo 

Rsh OW er = AL s effP N A J2 J dEP ( F*(E5 t )(<T BB (E> t ) 



' Eo 

^y^hc+NC+EK\^Vi i V) 



E h0 /E? 



dy 



) (17) 



where A is the detector area and L s e ^ = LD+0.3(km), as 
discussed below for shower rates, is the effective length of 
a detector of instrumental length L o . E%. is the primary 
neutrino energy, Eq has to be greater than or equal to 
the minimum energy at which the flux is known, and E^q 
is the minimum energy of the hadronic shower. The sum 
over 'i' is to account for different neutrino flavors. With 
the exception of v e CC interaction, for which we must 
set Eho = for the reason given below, we have to set 
Eho = Eo > ITeV as we do not know the flux below Eq 
and hence cannot account for all the showers produced 
below Eq. One needs E^q > ITeV for the showers to 
be detectable. For the black hole we assume that it has 
equal probability of decaying into any SM particles, and 
that it will always give a shower of energy around Eq and 
higher. 

Two important points: (i) We take the effective detec- 
tor length L s e ^ as the instrumental detector length Ljj 
•plus 0.3km. This is because, in addition to the showers 



produced inside the detecor, a conservative estimate is 
that shower signal produced 0.15km outside the detector, 
on any side, will easily reach the detector using a shower 
range of 0.3km for optical modules. This increases the 
shower rates by 30%. (ii) For showers from v e CC in- 
teraction, there is no need to set a lower limit on y as 
the electron energy will add to the hadron energy to con- 
tribute to the shower signal (e.g. effectively one can set 
Eho = 0)- This increases the shower rates dramatically 
— 30-50%. This is because in this case we can set min- 
imum y=0 and the CC cross section peaks around y=0 
giving a large percentage of the total shower events. 

(ii) Muon rates: CC interaction and v T CC in- 
teraction, with the tau decay r — > v T v^pL, both are the 
sources of muons |84|. For the former case, muon rate 
B£ c is given by 



R 



cc 



OO "n 

A P N A J dEl h F^{El h ) J dy- 

E 



^cc(EPy) 



L» ff (E^y)0( RiE ^ o) - {x^: ck +l sh ° w )\ 

(18) 

where E^q = Eq for the same reason as given above for 



iitrack 
r ■ 
mm 



muon track length required to detect a muon and l show 
~ 0.02fcm (<C x*^™ c ) is the typical shower size (for PeV 
energies) of the shower produced at the event vertex. 
The 9 function guarantees the exclusion of muons whose 
range is too small for them to be detected. L£^(2?JJ , y) 
is the effective detector length for muon rates: 



R(E^EP,y),E, l0 ) 



2x 



[itrack 



ishow 



where 



P PP 



■0E« 



■ fiE^Q 



(19) 



(20) 



is the average electromagnetic range, in a matter of den- 
sity p, of a muon of initial and final energies E^ and 
E^q, respectively Here a — 2.0MeVcm 2 /g accounts for 
the muon energy loss due to ionization and (3 — 4.2 x 
10~ 6 cTO 2 /g is due to pair production, bremmstrauhmg, 
and photonuclear energy losses 

We should emphasize that the effective detector length 
L^ff (E?^ , y) given above is the appropriate one for this 
case. Our definition of the effective detector length works 
for any value of the muon range while - R (- E ^- E ^°) works 

only for R ^ E "<>) > L D . 

Muon rate from the tau decay r — > v T v^pL is given by 
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E B^o/Bt 



Riz'Er.E^), 



L» ff {El T ,z')9 { ^Bil ) _ (a ££«* + ) , (21) 



where C T dec is defined in Eq^J E t0 = E for the 
same reason as given above for the shower and muon 
rate; dP d , ^ z •* gives the relevant decay distribution with 
E„ = z'E T = z'(l - y)ES T M 2000]; x^[ n ack and l show 
are the same as defined above. -k^// is defined in Eq. 
mil The integration over x gives total probability that 
a tau will decay with a tau decay length between and 



oo. The 6 function requires the muon range to be greater 
than the minimum distance required to detect a muon. 

(iii) Tau rates: We discuss two types of events that 
are unique to the presence of taus and v T 0> EE Il3 : 

(1) A tau produced in a v T CC interaction outside the 
detector decays (excluding the decay r — * Vr^iv^) inside 
the detector a track and a shower: 



OO V -y 

C f da" T (E p v) 

R T T dec (l shower) S 0.83,4pA^ y dE^F"* (E* t ) J dy 



En 

I L D +0.3(fem) 



-x/Cl ec {E T ) 



Vmin' 1C ' e + - 15 ( fem ) 



Lr)+0.3(fcm) 



; -x/£5 ec (B T ) 



(22) 



where 

Lleff( X ) = ( x ~ x min-0.15(fem)) 

L^// - (io + 0.15(fem) - <<™ cfe ), (23) 

where -E T o = Eq for the same reason as given above for 
the shower and muon rate. The 0.83 factor is to ex- 
clude the decay t — ► v^-paj^ which has a branching ra- 
tio of "0.17, and x r ^ ck ~ 0.25(fcm) is the minimum 
tau track length required to detect a tau. We have as- 
sumed that all the showers produced in these events will 
be detectable, so we do not need tau decay distribution 
function. This is a reasonable assumption as we will 
choose Eq = E t q — 0.5PeV which means almost all of 
the showers produced from tau decay will be above the 



detector threshold of ~ O.OOlPeV. The lower limit for 
x integration assures the tau decay length large enough 
for the tau to be separately detected from the shower. 
The expression given for these events in Ref . 01 includes 
some shower — track — shower events too. The expres- 
sion above gives only track — shower events by using 
the x dependent effective length where x is smaller than 
Lq + 0.3(fcm). Moreover, we have included the 0.3(km) 
in the x integration limit and 0.15(km) in the expres- 
sion for i^e//- These numbers follow from the reasoning 
given above in discussion of shower rates. 

(2) A tau produced in a v T CC interaction inside the 
detector decays (excluding the decay r — > v T fiv li ) inside 
the detector giving a shower, a track, and another shower. 
These are so called double bang events [7fij : 



1 _ SjO 

e?,_ L D +0.3(km) 

-x/C^(E T ) 



R^{2 shower) S 0.83ApN A J dE^F+fflJ J dy cc ^ ,y > J dx £r £g e// (a;),(24) 



E xZ$ r * ch +l" h ° 



where and all the other symbols have been defined above. 

L leff{x) = (Ld + 0.3(fcm) -x) (25) 
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TABLE III: Up and down events (yr^ 1 ) for E > 0.5PeV in 
the scenario 1:1:1; all upward events are integrated over nadir 
angle 8 < 84°; down showers are integrated over angle but 
muons and taus are not (see text for details). 





showers ( 


"P \ 
down ) 


muons ( 


_up_\ 

down ' 


taus (-j2E_) 

^ down / 




SM 


G2 


Gl 


SM 


G2 


Gl 


SM 


G2 


Gl 


WB 


2.6 
11 


2.7 
24 


3.1 

201 


3.0 
3.3 


2.7 

3.3 


i.i 

3.3 


0.11 
0.13 


.074 
0.13 


0.0086 
0.13 


SD 


163 


167 


202 


176 


165 


74 


4.7 


4.0 


0.54 


622 


74 S 


4725 


142 


142 


142 


4.6 


4.6 


4.6 


MB 


3.0 


3.1 


2.2 


6.4 


3.6 


0.66 


0.46 


0.20 


0.01 


30 


195 


1898 


18 


18 


18 


0.86 


0.86 


0.86 


PR 


32 


33 


34 


50 


38 


12 


2.4 


1.5 


0.13 


182 


534 


5284 


73 


73 


73 


3.5 


3.5 


3.5 



TABLE IV: Ratios of the ratios; here RRU 

^»^»2 taus down 



showers down 
■muons up 



and 







RR1 

RR1 


G2 




RR1 

RR1 


Gl 
SM 




RR:2 G2 


RH-2 G1 




1 : 2 


: 


1:1:1 


1 : 2 


: 


1 : 1 


: 1 






WB 


4.5 _ 
1.7 


2.7 




93 _ 
1.7 


56 


183 _ 
3.6 


=50 


1.5 


13 


SD 


2.1 
1.6 — 


1.3 




34 
1.6 


21 


64 
3.5 


18 


1.2 


8.7 


MB 


30 _ 
2.1 


14 


4.7 - LJ - 


1571 
2.1 


= 735 


2875 
4.7 


=610 


2.3 


46 


PR 


7.4 
1.6 — 


4.7 


3.6 _ 


249 
1.6 ~~ 


=157 


440 
3.6 - 


=121 


1.6 


18 



In addition to the above events for tagging taus, we 
have looked at the possibilty of detecting taus from the 
decay r — > v r \iv^ provided the tau decays inside the 
detector. These events have smaller rate than double 
bang or single shower events. Details of this calculation 
will be given elsewhere. 



VI. RESULTS AND DISCUSSION FOR EVENT 
RATES 

The results for event rates are summarized in the Ta- 
bles IHItlVll The event rates at angles below 60° turn out 
to be very small for current detectors, hence, we show 
events integrated to 9 = 84° nadir angle. 

In Table ITTT1 we give different upward and down events 
per year, in the scenario 1:1:1, for a lfcm 3 detector in ice. 
Our down shower events have been integrated over angle, 
however, the down tau and down muon events correspond 
to near horizon events. The showers contain events due 
to all neutrino flavors. These shower events in 1:1:1 sce- 
nario are 30- 40% larger than 1:2:0 scenario (not shown 
here). This is because, in the latter, we are excluding 
CC interaction and also the total flux is smaller due 
to the absence of tau regeneration effect. Upward muons 
(muon up) contain muons from CC interaction and the 
muonic tau decay. The muons in 1:1:1 scenario are about 
40% smaller than 1:2:0 scenario (not shown). They are 
not exactly 50% of the latter due to the contribution from 
tau decays. The upward taus (taus up) contain all three 
types of events described above (e.g. tau up = (track — 
shower) + (shower — track — shower) + r —* Vt^v^). 
The ratio (track — shower) / (shower — track — shower), 
not shown here, is very sensitive to the flux model and 
can be anywhere between 0.7 to 1.7. The r — > Vt-ixv^ 
tau events are always the smallest — less than 50% of 
the smaller of the single shower and double bang. While 
muons show differences between LSG (lTeV)and SM, the 
up-to-down showers and tau event ratios show a clearer 
differentiation between the two (Table ITTT^ . The number 
of events, though marginal for z/ T 's in WB and MB, are 
sufficient to make a clear distinction from SM for LSG 
(1 TeV) and distinction in some cases for LSG (2 TeV) 
in several years of data taking. The up tau events in 
the LSG (1 TeV) are especially severely suppressed as 



TABLE V: Same as Table|lTT|but with energy threshold 5PeV. 





showers ( 


up \ 
d.mim. J 


muons ( , up ) 


taus ( , up ) 

\ damn. / 




SM 


G2 


Gl 


SM 


G2 


Gl 


SM 


G2 


Gl 


WB 


0.26 


0.33 


0.15 


0.3 


0.14 


.002 


0.07 


0.04 


.001 


2.9 


13 


151 


0.84 


0.84 


0.84 


0.11 


0.11 


0.11 


SD 


11 


15 


8.7 


8.4 


5.6 


0.1 


2.5 


1.8 


0.05 


94 


188 


2714 


14 


14 


14 


3.1 


3.1 


3.1 


MB 


1.2 


1.4 


0.41 


2.0 


0.66 


.006 


0.41 


0.16 


.003 


21 


157 


1602 


9.2 


9.2 


9.2 


0.83 


0.83 


0.83 


PR 


6.5 


8 


2.9 


8.5 


3.7 


.04 


1.9 


0.93 


0.02 


82 


371 


4178 


26 


26 


26 


3.0 


3.0 


3.0 



compared to the SM. This again reflects the fact that 
tau decay is playing much weaker role in LSG (ITeV) 
(see Section IV). However, as expected, taus play a simi- 
lar role in SM and LSG(2TeV), though some suppression 
of v r is evident for LSG (2TeV) in Tables |Hl] and |IV| 
In Table ITVl we see the LSG/SM ratio of the ratios, as 
defined in the table caption, strongly differentiates SM 
from LSG(lTeV)in both of the flavor scenarios. In fact, 
these ratios have very weak dependence on the flavor sce- 
nario. They are spectacularly large in every flux model. 
Even LSG (2TeV) is clearly distinguished in all but the 
SDSS model. In the latter, the number of events is large 
enough that one may hope to discriminate between LSG 
(2TeV) and SM. 

We show similar tables with an energy threshold at 
5PeV (TableslVland IVI|) . Here we see the tau events did 
not change much in SM and LSG (2 TeV) model. How- 
ever, in LSG (lTeV)up tau events have decreased by an 
order of magnitude, generally with enough down events 
to make them even more useful in differentiating between 
the two at 5PeV as compared to 0.5PeV. Although not 
shown here, by comparing 5PeV and 0.5PeV results for 
showers and muons one expects that LSG(2 TeV) model 
around lOPeV thresholds will differ from SM to the ex- 
tent that LSG (1 TeV) does from SM around 0.5PeV. 
However, the event rates may be too small to do a sta- 
tistical analysis. 

Looking at the event rates for different flux models 
in Table InTl we see the PR and SDSS flux models pro- 
vide large enough events to do a statistical analysis. For 
WB and MB models, only the down shower rates in 
LSG(lTeV) are large enough to differentiate it from the 
SM by looking at the up-to-down shower ratios. 
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TABLE VI: Same as Table |TV| but with energy threshold 
5PeV. 









G2 


HM 


Gl 


HH l 2 G 2 


HH'2 G i 






nm 




RR1 


'M 


-R-R2.9 M 


RR2s M 




1 : 2 


: 


1:1:1 


1:2:0 


1:1:1 






WB 


11 




9.3 


8400 


7826 


1.7 


74 


SD 


3.2 




3.0 


2481 


1419 


1.4 


48 


MB 


27 




23 


3.0e4 


2.7e4 


2.6 


161 


PR 


12 




10 


1.2e4 


1.0e4 


2.0 


100 



VII. SUMMARY OF RESULTS AND 
CONCLUSION 

We found complete numerical solutions to the system 
of coupled equations that include the most important ef- 
fects for transport through Earth of v e , v^, v r and t fluxes 
above 0.5 PeV. In Fig. 4, we presented results of angular 
distributions of total neutrino flux in our example of the 
diffuse flux model by R. J. Protheroe ^l| 1 however, the 
qualitative features of this figure are common to the other 
flux models (K. Mannheim (B) \v^, Waxman Bahcall 
[To|. and SDSS 0). Fluxes in this figure are integrated 
from 0.5 PeV upward, showing the effects of including 
low scale gravity enhancement to the lepton deep inelas- 
tic cross sections, with no v T and full v T mixing into the 
incident flux. This figure also show that the v T regen- 
eration from t decay enhances the "through earth" , or 
"upward" fluxes significantly more in the standard model 
than in the models with low scale gravity enhancements 
included, as seen on the curves where v T is mixed into 
the flux incident on Earth. The standard model flux is 
obviously higher at nadir angles smaller than 80 degrees, 
while the differences between the fluxes with standard 
model interactions only and those with low scale gravity 
included are much larger at small nadir angles in the case 
that v T is mixed into the incident flux than in the case 
when there it is not. Next in Figs. 5-8, we showed the 
equivalent angular distributions for the v T flux alone to 
emphasize the observation just summarized, that is, as 
compared to v e and v T can serve better to differenti- 
ate between SM and LSG at energies below lOPeV. 

As established by the angular distribution graphs, the 
qualitative features are shared by all the models, so we 
gave only the Protheroe model results in plotting the 
energy distribution of flux integrated over angles in the 
range from 0.5 PeV to 20 PeV in Figs. 9 and 10. These 
plots show that in this energy range, the low scale gravity 
interactions rapidly suppress the upward flux compared 
to the standard model. They also indicate the fact that 
the regeneration of v T flux is much less significant when 
low scale gravity is turned on, as clearly indicated by the 
flux difference curves in the energy range between 0.5 
PeV and 2.0 PeV. 

In the series of graphs from Figs. 11 through 14, we 
displayed the three-dimensional plots of the total and 
v T — only fluxes for the standard model and for the low 
scale gravity, M = 1 TeV case. These indicate in de- 
tail where the maximum flux differences are in angle and 



energy. 

Next we looked at the flux ratios Rl, R2, and R3 as 
defined in Eqs. ED The results are given in Tables H| and 
ITT1 Here again the distinction between v T and v^, + v e 
fluxes is evident. The distinctions among flux models 
are largely washed out by the LSG dynamics at higher 
energies. For example, in Table [HJ Rl and R2 are the 
same for all the flux models given here. 

In Sections V and VI we presented the defining equa- 
tions for our shower, muon and tau rates and the results 
of our rate calculation. The story is summarized in Ta- 
bles III- VI. Using a cutoff of 0.5PeV, we found that the 
events rates in showers and muon categories are large 
enough to make meaningful statements about the dis- 
tinction between SM and LSG with a 1 TeV, in all flux 
models and both flavor scenarios with a 2-3 years of run- 
ning. An interesting feature of the LSG (2TeV) entries in 
Tables ITTT1 and HVl is that the down shower events may be 
enhanced enough compared to SM to distinguish between 
the two in WB, MB, and PR, and possibly SD too. The 
ratio of ratios in Tables IIVI and IVII compares the LSG 
shower down/muon up ratios to the ones for SM. Table 
IV and VI also show the same for taus down/taus up. 
This diagnostic is especially sensitive to the difference 
between SM and LSG. It also shows us that this ratio 
of the ratios for showers and muons is almost the same 
in both of the flavor scenarios. The importance of the 
taus in differentiating between SM and LSG (1 TeV) is 
realized by looking at the tagged tau events (Tables IV 
and VI). For tagged tau events, the difference between 
LSG (1 TeV) and the SM varies from an order of mag- 
nitude to two orders of magnitude for energy thresholds 
of 0.5-5 PeV. However, we caution the reader again that 
the statistics are low in this case. Basically the tau story 
can be summarized by saying that any upward tau event 
establishes (1) v T 's presence in the neutrino flux incident 
on earth (2) exclusion of LSG with 1 TeV scale or any 
model of enhanced cross section of comparable size in the 
1-10 PeV range. 

The 5 PeV threshold results in Table V and VI show 
the same patterns as in the 0.5PeV tables. The distinc- 
tion between SM and LSG (2 TeV) are now sharper in 
the ratios, but the statistics in some cases are low, so 
that one needs to have 5-10 years of data to draw strong 
conclusions that apply to all flux models. 

We conclude on the basis of our flux and event rate 
study that with a threshold of 0.5 PeV, the shower and 
muon event ratios have sufficient events in all flux and 
lepton flavor models to make clear distinctions between 
SM and LSG with a mass scale 2 TeV and below. Going 
above 2 TeV, one finds that whether distinctions can be 
made depende upon the flux model. The situation is not 
so clear. Because the requirements on v r identification 
are so stringent, only a few events to a fraction of an event 
will be expected, depending upon flux model, up or down 
event and LSG scale value. One point is perfectly clear: 
any upward tau event excludes LSG with a scale around 
ITeV. 



14 



Given the intense experimental activity in the field, 
we expect that data will yield many insights in the com- 
ing decade when analysed with techniques like the ones 
presented here. 
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